MIT  LIBRARIES  DUPL 


3  9080  02245   1980 


28 

14  M 

1 


DEWEY 


The  Impact  of  novel  Treatments  on  A6  Burden  in  Alzheimer's 
Disease:  Insights  from  a  Mathematical  Model 

David  L.  Craft,  Lawrence  M.  Wein,  Dennis  J.  Selkoe 

Sloan  Working  Paper  Number  4 1 56 
January  31,2001 


3 


The  Impact  of  Novel  Treatments  on  A/5  Burden 
in  Alzheimer's  Disease:  Insights  from  a  Mathematical  Model 

David  L.  Ciaft\  Lawrence  M.  Wein"^,  Dennis  J.  Selkoe^* 

^Operations  Research  Center,  MIT.  Cambridge,  K4A  02139 

'^ Sloan  School  of  Management,  MIT,  Cambridge,  MA  02139 

Center  for  Neurologic  Diseases,  Harvard  Medical  School,  Brigham  and  Women's  Hospital, 

Boston,  MA  02115 


'Person  to  whom  correspondence  should  be  addressed:  email:  si'lkoe fend. bwh. harvard. eihi 


Abstract 

Motivated  by  recent  therapeutic  initiatives  for  Alzheimer's  disease,  we  developed  a 
mathematical  model  of  the  accumulation  of  amyloid  /3-protein  (A/i)  in  the  brain.  The 
model  incorporates  the  production  and  clearance  of  kl5  monomers,  and  the  elongation 
and  fragmentation  of  polymers  by  monomer  aggregation  and  break-off,  respectively. 
Our  analysis  suggests  that  A/3  dynamics  are  dictated  by  a  single  tmitless  measure  re- 
ferred to  as  the  polymerization  ratio,  which  is  the  product  of  the  production  and  elon- 
gation rates  divided  by  the  product  of  the  clearance  and  fragmentation  rates.  Cerebral 
k(i  burden  attains  a  finite  steady-state  level  if  this  ratio  is  less  than  one,  and  undergoes 
sustained  growth  if  this  ratio  is  greater  than  one.  The  highly  nonlinear  relationship  be- 
tween the  polymerization  ratio  and  the  steady-state  A/3  burden  implies  that  a  modest 
reduction  in  the  polymerization  ratio  (e.g.,  a  40%  reduction  in  the  production  rate  by 
a  7-secretase  inhibitor)  achieves  a  significant  (e.g.,  18-fold)  decrease  in  the  A/3  burden. 
Our  model  also  predicts  that  after  initiation  or  discontinuation  of  treatment,  it  may 
take  months  to  reach  a  new  steady-state  A/3  burden.  Taken  together,  our  findings 
suggest  that  the  research  comnumity  should  focus  on  developing  agents  that  provide 
a  modest  reduction  of  the  polymerization  ratio  while  avoiding  long-term  toxicity.  Fi- 
nally, our  model  can  be  used  to  indirectly  estimate  several  crucial  parameters  that 
are  difficult  to  measure  directly:  the  production  rate,  the  fragmentation  rate  and  the 
strength  (i.e.,  percentage  inhibition  or  enhancement)  of  treatment. 


Introduction 

The  cerebral  build-up  of  amyloid  /^-protein  (A/3)  and  its  aggregation  into  oligomers  and 
polymers  are  key  components  of  the  pathogenesis  of  Alzheimer's  disease  (AD)  '.  Recent 
studies  suggest  that  A/3  accumulation  plays  an  important  role  in  AD  neurodegeneration 
■^'■^,  and  a  variety  of  A(5  treatment  strategies  are  being  aggressively  pursued,  including  the 
enhancement  of  A/?  clearance  via  immunization  with  A/3  ^'^^  and  the  reduction  of  A/3 
production  by  inhibiting  7-secretase,  which  is  one  of  the  two  enzymes  that  cuts  the  amyloid 
precursor  protein  (APP)  to  release  the  A/3  fragment  ^•^.  The  kinetics  of  A/3  production, 
aggregation  (polymerization),  disaggregation  and  clearance,  and  the  effect  of  these  novel 
treatments  on  the  time-dependent  variation  of  A/3  levels  are  not  well  understood.  Although 
there  exist  several  mathematical  models  that  focus  on  either  fibrillogenesis  ^"^^  or  plaque 
formation  ^''^'^^  these  studies  neither  consider  the  impact  of  treatment  nor  allow  for  a 
contiiniously  renewable  source  of  A/3  from  APP  molecules.  Here  we  formulate  and  analyze  a 
simple  mathematical  model  that  tracks  the  dynamics  of  A/3  production,  polymer  elongation, 
fragmentation  and  clearance  during  the  course  of  treatment.  Simple  formulas  and  numerical 
results  are  presented  that  elucidate  the  impact  of  treatment  on  A/3  burden. 

Mathematical  model 

In  AD  brain  tissue,  the  level  of  extracellular  (and  perhaps  intracellular)  k[5  (particularly 
A/342,  the  42  amino-acid  form  of  A/3)  increases  over  time,  which  gives  rise  to  polymeriza- 
tion. Some  of  these  polymers  cluster  into  light  microscopically-visible  particles  consisting 
of  A/342,  which  can  accumulate  to  create  diffuse  (amorphous)  platiues,  the  apparent  initial 
neuropathological  lesion  of  AD.  Further,  some  polymers  containing  A/^42  and/or  A/340  t'^n 
eventually  fold  into  long  hlamentous  assemblies  called  amyloid  hbrils,  and  clumped  masses 
of  these  are  referred  to  as  amyloid  (senile)  plaques.  Our  mathematical  model  idealizes  this 


process  in  several  iinjxjrtaiit  ways.  First,  we  only  consider  extracellular  A/?,  which  is  likely  to 
be  in  a  dynamic  equilibrium  with  intracellular  A/i.  We  also  do  not  distinfi,uish  between  A/i4o 
and  A/?42,  or  between  soluble  and  insoluble  AS,  because  there  is  insufficient  information  in 
the  current  literature  about  the  precise  dynamic  changes  in  these  species  over  time  in  the 
preclinical  and  clinical  phases  of  AD.  Finally,  because  our  primary  goal  is  to  understand 
the  impact  of  treatment  on  total  cerebral  A/3  burden  (as  opposed  to  the  number  and  size 
of  plaques),  we  restrict  our  attention  to  A/3  polymerization,  and  do  not  attempt  to  capture 
the  downstream  processes  of  fibrillization  or  plaque  formation  in  vivo.  However,  later  in 
the  paper  we  discuss  how  this  model  might  be  interpreted  in  terms  of  -  and  generalized  to 
incorporate      plaque  formation. 

Our  model  is  an  infinite  system  of  nonlinear  differential  equations  that  tracks  the 
temporal  evolution  of  the  concentration  of  extracellular  A/3  i-mers,  for  2  =  1,2,....  It  is  a 
variant  of  the  Smoluchowski  equation  -'^^~\  which  has  been  used  for  nearly  a  century  to  study 
aggregation  processes  such  as  polymerization,  galaxy  formation,  crystallization,  and  cloud 
formation  2-2-23  vVfhiie  this  nrodel  typically  considers  aggregation  in  closed  systems,  the 
novelty  of  our  model  is  to  simultaneously  incorporate  fragmentation,  a  source  (A/i  monomer 
production)  and  a  sink  (A/i  monomer  loss). 

For  t  =  1.  2, . . . ,  we  let  c,(f)  denote  the  concentration  of  A/?  z-mer  at  time  t\  i.e.,  i  =  1 
denotes  monomer,  t  =  2  denotes  dimer,  etc.  The  differential  equations  dictate  the  time  rate 
of  change  of  r,(i),  denoted  by  Cj(t),  and  are  given  by 
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The  model  captures  four  processes  -  production,  elongation,  fragmentation,  and  loss 
(i.e.,  clearance)  -  and  for  ease  of  reference  the  corresponding  mnemonic  parameters  are 
displayed  in  Table  1.  In  equation  (1),  A/3  monomers  are  produced  via  cleavage  of  the  APP 
at  a  constant  rate  p.  and  are  lost  at  time  t  at  rate  lci{t).  That  is,  A/3  monomers  live  for  /"^ 
time  units  on  average  before  being  cleared  by  cell  internalization  (e.g.,  microglial  ingestion), 
degradation  by  proteases  or  removal  from  the  brain  via  the  circulation. 


Parameter 

Description 

Value 

P 

Production  rate  of  A/3  monomers 

1.02  X  10-1^  M  sec-i 

I 

Loss  rate  constant  for  monomers 

1  hr^^ 

e 

Elongation  rate  constant  (only  monomer  addition) 

90  M-'  sec"' 

f 

Fragmentation  rate  constant  (only  monomer  break-off) 

3.93  X  10"^  sec-' 

Table  1:  Some  theoretical  parameter  values  for  the  model,  which  yield  a  polymerization  ratio 

of  r  =  pe/(fl)  =  0.84. 

Equations  (l)-(3)  assume  that  elongation  occurs  by  monomer  addition  with  elongation 
constant  e.  so  that  eci(f)Cj_i(t)  is  the  rate  at  which  (?  —  l)-mers  elongate  to  z-mers  for  ;  >  1. 
Similarly,  to  limit  the  number  of  model  parameters  recjuired,  we  assume  that  fragmentation 
occurs  only  by  monomer  break-offs,  with  fragmentation  rate  constant  /,  so  the  fci{t)  is  the 
rate  at  which  an  ;  — mer  fragments  into  an  {i  —  l)-mer  and  a  monomer.  Hence,  for  i  >  2,  the 
concentrations  of  /'— mers  in  equations  (2)-(3)  increase  due  to  elongation  of  (i  -  l)-mers  and 
fragmentation  of  (/  +  l)-mers,  and  decrease  due  to  elongation  and  fragmentation  of  /-mers. 


The  factor  2  in  equation  (1)  arises  because  the  elongation  to  a  dimer  requires  two  monomers, 
and  the  term  //2  in  equation  (2)  stems  from  the  fact  that  a  dimer  can  only  fragment  in 
one  location,  while  larger  «-mers  possess  two  potential  fragmentation  sites.  Fragmentation 
of  dimers  occurs  at  rate  //2  and  creates  two  monomers,  and  hence  the  factors  2  and  1/2 
cancel  each  other  out  in  the  fco  term  in  ecjuation  (1). 

An  alternate  modeling  approach  is  to  allow  direct  clearance  of  monomers  off  of  i-mers 
(e.g.,  to  represent  microglial  ingestion  of  polymers),  rather  than  requiring  a  two-step  proce- 
dure of  monomer  fragmentation  followed  by  monomer  clearance.  This  alternative  would  re- 
sult in  the  omission  of  the  "all  fragmentations"  term  in  equation  (1)  and  the  re-interpretation 
of  /  as  an  ingestion  rate,  but  would  not  change  the  qualitative  nature  of  our  results:  it 
would  mainly  change  the  precise  conditions  for  a  steady-state  solution  and  the  relative  val- 
ues of  steady-state  monomer  and  dimer  concentrations,  and  slightly  hasten  the  approach  to 
a  post-treatment  steady-state.  Finally,  while  AjS  treatment  parameters  are  not  explicitly 
incorporated  into  equations  (l)-(3),  treatment  can  be  included  in  the  model,  as  explained 
later. 

Results 

Steady-state  solution 

Our  analysis  begins  by  seeking  a  steady-state  (i.e.,  equilibrium)  solution,  Cj,  to  equations  (1)- 
(3).  To  describe  our  results  in  a  concise  manner,  we  define  the  polmerization  ratio 

which  is  a  unitless  quantity  that  incorporates  the  four  model  parameters.  Setting  the  left 
side  of  equations  (l)-(3)  (i.e.,  the  rates  of  change  of  i-mer  concentration)  to  zero  and  solving 
for  Ci  reveals  that  there  are  two  regimes:  a  steady-state  (or  subcritical)  regime  where  /■  <  1 
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and  a  supercritical  regime  where  r  >  1.  In  the  former  case,  the  A/?  levels  settle  into  a 
steady-state  where 

d  =  y,      and     c,  =  2^r'-'     for     /  >  2.  (5) 

Hence,  ?— mer  concentrations  decay  geometrically  for  i  >  2,  and  monomers  are  more  preva- 
lent than  dimers  if  r  <  0.5,  and  less  prevalent  than  dimers  if  r  >  0.5.  In  contrast,  if  r  >  1 
then  the  A/j  burden  grows  indefinitely.  Because  the  A/?  burden  appears  to  remain  relatively 
stable  over  time  in  symptomatic  AD  patients  ^-^-'^s,  we  focus  primarily  on  the  steady-state 
regime,  and  defer  until  later  a  discussion  of  the  supercritical  regime. 

Equation  (5)  implies  that  the  steady-state  value  for  the  total  A/3  concentration  (i.e., 
total  number  of  Af3  molecules  in  the  system,  whether  they  exist  as  monomers  or  ?'-mers). 
which  is  denoted  by  r  =  ^j^i  ic^  and  referred  to  as  the  Ap  burden,  is  given  by 

Notice  that  the  total  Aft  burden  c  approaches  infinity  as  the  polymerization  ratio  r  appoaches 
1.  We  now  use  equations  (5)-(6)  to  estimate  the  values  of  the  model  parameters,  and  these 
parameter  values,  in  turn,  allow  us  to  calibrate  the  key  relationship  (6)  (i.e.,  A/i  burden)  to 
the  clinical  setting. 

Parameter  estimation 

The  value  of  the  A;:^  monomer  loss  rate  /  in  Table  1  corresponds  to  a  half-life  of  41.6  minutes, 
which  is  close  to  the  crude  experimental  value  of  38  minutes  reported  by  one  group  in  the 
brains  of  APP  transgenic  mice  ^.  Protofibrils  ~',  fibrils  ''■  "  and  plaques  ^^  appear  to  grow 
primaril\-  via  A.J  monomer  addition,  and  we  use  an  elongation  rate  e  in  Table  1  taken  from 
a  synthetic  fibril  analysis  ^'*,  which  is  within  a  factor  of  two  of  other  estimates  for  fibrils  ^^ 
and  protofibrils  ^' . 
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The  fragmentation  rate  /  and  the  monomer  production  rate  p  are  difficult  to  estimate 
empirically.  Fortunately,  we  can  use  equations  (5)-{6),  together  witli  A/i  estimates  from  hu- 
man AD  brains,  to  estimate  these  two  quantities  as  follows.  A  rec(Mit  study  *  estimated  that 
1.4%  of  A/?  in  the  human  brain  is  soluble,  and  found  that  this  soluble  compartment  consists 
primarily  of  monomers,  dimers  and  trimers.  A  related  study  '-^^  measured  the  weighted  av- 
erage (of  A/i4o  and  A/?42)  soluble  fraction  of  A/i  in  human  brains  to  be  1.2%.  By  taking  the 
average  of  these  two  estimates  (we  ignore  a  third  estimate  that  is  only  0.05%  '^'')  and  assuming 
that  the  soluble  dimers  and  trimers  are  exactly  offset  by  any  insoluble  monomers,  we  equate 
the  fraction  of  A/?  that  consists  of  monomers,  ci/c,  to  0.013.  We  recognize  that  this  analysis 
will  underestimate  the  percentage  of  soluble  A/3  in  the  brain  if  there  are  large  amounts  of 
soluble  dimers  and  trimers  in  AD  brain.  By  equations  (5)-(6),  Ci/c  =  (2/(1  -  r)''  -  1)^\ 
and  setting  this  expression  equal  to  0.013  yields  the  polymerization  ratio  r  =  0.84.  We  also 
equate  the  A/i  burden  on  the  right  side  of  equation  (6)  to  the  average  of  the  total  A/3  levels 
(A/3.}o  plus  A/342)  found  in  five  different  cortical  regions  of  wet  brain  tissue  of  patients  with  a 
clinical  dementia  rating  (CDR)  score  of  5.0  (severe  dementia)  in  Table  1  of  -  (this  estimate 
of  2819  pmol/g  is  similar  in  magnitude  to  those  in  ■'^),  by  assuming  that  the  density  of 
wet  cortical  tissue  (which  is  about  87%  water)  is  equal  to  that  of  water.  Substituting  our 
estimates  for  /  and  r  into  the  right  side  of  equation  (6)  and  equating  this  expression  to  2819 
pmol/g  yields  the  value  of  the  production  rate  p  in  Table  1.  Finally,  because  r  =  pe/{fl), 
the  polymerization  ratio  estimate  r  =  0.84,  together  with  our  earlier  estimates  of  e,  I  and  p, 
yield  the  value  for  the  fragmentation  rate  /  in  Table  1.  Note  that  this  value  of  /  is  about  70 
times  smaller  than  the  value  of  the  loss  rate  /,  which  suggests  that  within  the  context  of  our 
model  an  A/3  vaccine  ^'  ^  that  directly  ingests  monomers  off  of  oligomers  can  be  accurately 
represented  as  a  fragmentation  enhancer. 


A/3  burden 

With  these  parameter  values  in  hand,  we  now  return  to  equation  (6)  for  the  steady-state  total 
A/3  burden  c.  This  nonlinear  relationship  between  the  polymerization  ratio  and  the  steady- 
state  A/?  burden  is  shown  in  Figure  1,  and  implies  that  a  given  treatment  will  generally 
achieve  a  greater  (relative  and  absolute)  reduction  in  steady-state  A/?  burden  for  patients 
with  higher  pretreatment  A/3  burdens.  To  position  this  relationship  within  the  clinical 
context,  we  convert  the  A/3  burden  in  Figure  1  to  the  CDR  scale  using  the  data  in  Table  1  of 
^.  This  conversion  allows  us  to  explore  the  potential  clinical  impact  (i.e.,  assuming  that  the 
neurodegeneration  associated  with  A/9  can  be  reversed)  of  a  reduction  in  the  polymerization 
ratio  via  treatment.  Recent  studies  reveal  that  A/?  peptide  vaccination  reduces  not  only 
A/3  burden  but  also  cognitive  impairment  in  APP  transgenic  mice  ^~-  ",  implying  that  the 
neurodegeneration  is  at  least  partially  reversible  in  mice.  We  note  that  the  CDR  scales  on 
the  right  ordinates  of  Figures  1,  2  and  3  can  be  ignored  if  the  reader  believes  that  CDR  will 
not  fall  in  synchrony  with  A/3  burden  in  humans;  the  latter  ordinate  is  still  valid. 

7-secretase  inhibitors 

Because  of  their  potentially  adverse  effect  on  Notch  signalling  34-3(^  7-secretase  inhibitors 
might  best  be  used  to  reduce  only  modestly  the  A/^  monomer  production  rate,  perhaps 
by  20-6U%.  To  assess  the  effect  of  such  treatment,  we  numerically  solve  equations  (l)-(3) 
assuming  that  a  symptomatic  patient  is  in  a  pretreatment  steady-state  on  days  0-250,  and  is 
administered  a  7-secretase  inhibitor  that  decreases  the  A/j  production  rate  p  by  40%  on  days 
250-620  (i.e.,  for  1  year).  As  shown  in  Figure  2,  the  post-treatment  steady-state  represents 
an  18-fold  drop  in  the  A,3  burden  and  a  theoretical  decrease  in  the  concommitant  CDR  score 
from  5.0  (severe  dementia)  to  0.0  (no  dementia).  The  perturbation  dynamics  in  Figure  2 
are  rather  sluggish:   the  relaxation  time,  which  we  define  to  be  the  time  from  the  start  of 
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treatment  until  the  A/?  burden  drops  to  90%  of  the  way  towards  the  post-treatment  steady- 
state  (i.e.,  tlie  time  until  the  A/i  burden  equals  the  steady-state  post-treatment  A/?  burden 
plus  0.1  times  the  difference  of  the  steady-state  pretreatment  and  post-treatment  burdens) 
is  113  days.  Interestingly,  it  takes  much  longer  ~  about  600  days  after  the  discontinuation  of 
treatment  for  the  A/3  burden  to  revert  to  90%  of  its  pretreatment  steady-state  (assuming 
a  polymerization  ratio  r  of  0.84  (Table  1)). 

Comparing  different  treatments 

Equation  (6)  allows  us  to  provide  a  comparison  across  various  treatment  approaches.  Witliin 

the  context  of  our  model,  treatments  can  affect  all  four  parameters:  7-secretase  inhibitors 

^'  **,  or  other  agents  that  target  the  upstream  mechanism  by  which  A/3  is  produced  from 

APP  (e.g.,  /j-secretase  inhibitors),  reduce  the  production  rate  p;  agents  that  promote  the 

fragmentation  of  monomers  from  polymers  and/or  enhance  the  monomer  clearance  rate  ■*•  ^ 

increase  the  fragmentation  rate  /  and  loss  rate  /,  respectively;  and  agents  that  inhibit  the 

deposition  of  A/3  monomers  reduce  the  elongation  rate  e.   By  equations  (4)  and  (6),  if  the 

production  rate  p  is  reduced  by  x%,  then  the  resulting  A/^  burden  c  is  x%  lower  than  if 

the  elongation  rate  e  is  reduced  by  x%;.   Similarly,  if  the  fragmentation  rate  /  is  increased 

by  y%,  then  the  resulting  A/3  burden  is  y%  higher  than  if  the  loss  rate  is  increased  by  y%. 

Hence,  the  parameters  p  and  /  have  somewhat  more  leverage  than  e  and  /,  respectively, 

in  reducing  the  A/3  burden.  Moreover,  an  x%  reduction  in  the  production  rate  (elongation 

rate,  respectively)  has  the  same  effect  as  a  y%  increase  in  the  loss  rate  (fragmentation  rate, 

respectively),  where 

lOOx 
^~  100 -x" 

For  example.  Figure  2,  which  was  computed  using  a  7-secretase  inhibitor  that  reduces  the 

production  rate  by  40%i,  would  also  resvilt  from  a  treatment  that  increased  the  loss  rate  by 
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66.7%. 

Unfortunatt'ly,  there  is  no  simple  formula  to  compare  the  A/?  burden  reduction  achieved 
by  changes  in  p  vs.  /  or  by  changes  in  e  vs.  /.  However,  in  Figure  3  we  plot  th(>  post- 
treatment  steady-state  Af3  burden  (and  corresponding  CDR  score)  as  a  function  of  Ijoth 
the  percentage  reduction  in  production  rate  p  achieved  by  a  7-secretase  inhibitor  and  the 
percentage  increase  in  fragmentation  rate  /  achieved  by,  e.g.,  an  A/i  vaccine  (recall  that  /  is 
70  times  larger  than  /  in  Table  1).  This  graph  strongly  suggests  tliat  only  a  modest  (e.g., 
10%)  change  in  these  parameter  values  is  recjuired  to  achieve  a  several-fold  reduction  in  A(3 
burden  and  potentially  a  clinically  significant  effect  on  dementia.  Because  A/3  monomers 
need  to  be  cleared  after  they  break  off,  this  figure  also  shows  that  an  infinite  fragmentation 
rate  still  leads  to  a  finite  steady-state  A/3  burden,  given  the  limits  of  the  monomer  loss  rate 
/. 

Additional  computational  results  (not  shown  here)  reveal  that  a  fragmentation  rate 
enhancer  has  a  smaller  (i.e.,  faster)  relaxation  time  than  a  production  rate  inhibitor,  probably 
because  the  relaxation  time  depends  in  large  part  on  the  fragmentation  rate.  For  example, 
an  agent  that  increases  the  fragmentation  rate  by  118%!  achieves  the  same  post-treatment 
steady-state  A/3  burden  as  the  40%  7-secretase  inhibitor  in  Figure  2,  but  the  relaxation  time 
is  only  40  days  compared  to  113  days  for  the  7-secretase  inhibitor. 

Supercritical  regime 

Thus  far,  we  have  assumed  that  /•  <  1.  so  that  a  steady-state  is  achieved.  If  r  >  1  then 
the  A/3  burden  grows  indefinitely,  eventually  (perhaps  after  many  years)  increasing  linearly 
at  rate  p{r  —  l)/r.  Moreover,  for  all  r  >  1,  the  distribution  of  polymers  no  longer  decays 
geometrically,  but  tends  over  time  to  a  uniform  distribution,  where  each  ?-mer  (starting  with 
smaller  values  of  t)  successively  achieves  the  concentration  ci  =  //e,  c,  =  2//e  for  i.  >  2. 
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Heme,  there  are  three  possibilities  for  the  effect  of  treatment:  (i)  the  A/3  burden 
is  in  a  pretreatment  steady-state  and  drops  to  a  lower  post-treatment  steady-state,  as  in 
Figure  2;  (ii)  treatment  causes  the  A/i  burden  to  shift  from  the  supercritical  regime  to  the 
steady-state  regime  (Figure  4a);  and  (iii)  treatment  does  not  allow  the  A.i  burden  to  exit 
the  supercritical  regime,  although  it  does  lower  the  growth  rate  (Figure  4b).  Comparison 
of  Figure  4a  to  Figure  2  shows  that  it  may  take  nmch  longer  to  achieve  a  post-treatment 
steady-state  in  case  (ii)  than  in  case  (i). 

Sensitivity  analysis 

Because  we  do  not  possess  precise  estimates  for  our  parameter  values,  we  performed  a  sen- 
sitivity analysis  that  varies  the  pretreatment  polymerization  ratio  /•,  which  is  the  primary 
driver  of  the  A/3  dynamics,  according  to  our  analysis.  We  use  the  same  parameter  estimation 
technique  as  before:  fix  /  and  e  to  their  values  in  Table  1,  find  r  from  the  fraction  of  A/? 
monomer  ci/c,  use  the  average  A/?  burden  of  2819  pmol/g  '^  to  derive  p,  and  use  r  and  p  to 
calculate  /.  Rather  than  eciuating  Ci/c  to  1.3%,  we  now  allow  it  to  vary  in  order  to  generate 
a  wide  range  of  values  for  r,  which  in  turn  causes  p  and  /  to  change.  In  Figures  5a  and  5b, 
we  investigate  the  effect  caused  by  a  40%  reduction  in  the  production  rate  (via  a  7-secretase 
inhibitor)  as  a  fvmction  of  the  pretreatment  polymerization  ratio  7".  Figure  5a  plots  the  pre- 
treatment steady-state  A/3  burden  divided  by  the  post-treatment  steady-state  A/3  burden 
(i.e.,  the  fold-reduction  in  A/3  burden),  and  Figure  5b  shows  the  relaxation  time,  which  was 
defined  earlier.  Recalling  our  assumption  that  a  three-fold  reduction  in  AI3  burden  leads  to 
a  clinically  significant  improvement  in  CDR  score  (Figure  1),  Figure  5a  shows  that  the  40% 
production  inhibition  appears  capable  of  a  significant  clinical  improvement  for  essentially  all 
practical  values  of  r  (i.e..  corresponding  to  percentages  of  monoiueric  A/3,  ci/c,  less  than 
30%),  and  a  l-to-2  log  drop  in  A,3  burden  as  r  increases  from  0.76  to  0.93.  Figure  5b  shows 
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that  our  base  case  of  r  =  0.84  is  on  the  steep  portion  of  the  curve  of  relaxation  time  versus 
polymerization  ratio.  Hence,  if  the  fraction  of  monomeric  A/3  is  actually  5%  rather  than  1% 
of  total  A/i,  then  r  drops  from  0.84  to  0.7  and  the  relaxation  time  decreases  from  113  days 
to  only  two  weeks. 

Discussion 

Despite  the  existence  of  some  mathematical  models  for  A/i  fibrilli/ation  and  plaque  growth, 
this  paper  appears  to  represent  the  first  attempt  to  use  a  mathematical  model  to  assess 
the  effect  of  treatment  on  A/3  burden  in  Alzheimer's  disease.  Our  model  is  purposefully 
parsimonious,  and  contains  only  four  parameters:  A/3  monomer  production  rate  (p),  A/3 
monomer  loss  rate  (/),  polymer  elongation  rate  (e),  and  polymer  fragmentation  rate  (/). 

The  most  basic  result  of  our  analysis  is  that  there  are  two  possible  regimes,  depending 
upon  the  value  of  the  polymerization  ratio  r,  which  equals  pe/{fl).  If  this  unitless  quantity  is 
less  than  1,  then  the  A/9  burden  enters  a  steady-state  regime  and  takes  on  a  finite  value.  If  the 
polymerization  ratio  is  greater  than  one,  then  the  A/3  burden  grows  indehnitely,  eventually 
increasing  linearly  at  rate  p{r  -  l)/r.  Researchers  in  other  disciplines  have  found  that 
the  dynamics  of  their  complex  systems  are  dictated  by  a  unitless  measure  that  leads  to  a 
subcritical  regime  if  the  measure  is  less  than  1  and  a  supercritical  regime  if  the  measure  is 
greater  than  1;  two  examples  are  the  basic  reproductive  rate  Rq  of  an  infectious  pathogen 
in  epidemiology  ^^  and  the  traffic  intensity  p  in  queueing  theory  ^^.  Moreover,  as  with  AD, 
in  these  disciplines  the  ultimate  measures  of  performance  (size  of  the  epidemic  and  queue 
length,  respectively)  have  a  highly  nonlinear  relationship  with  the  unitless  measure.  Applied 
researchers  in  these  disciplines  have  adopted  these  unitless  measures  as  a  central  part  of 
\hc\r  mental  model  of  these  complex  systems,  and  routinely  estimate  them  and  monitor 
achieved  reductions  in  them.   We  propose  that  AD  researchers  follow  suit  by  adopting  the 
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polymerization  ratio  pe/{fl)  as  the  key  "A/3  driver"  in  this  disease. 

Several  studies  have  shown  that  the  A/3  burden  does  not  appear  to  be  tightly  correlated 
with  th(>  duration  or  clinical  severity  of  AD  21-26^  which  suggests  that  the  polymerization 
ratio  is  less  than  1  for  symptomatic  AD  patients.  Although  previous  mathematical  models 
of  plaque  formation  '*'  '^  have  also  proposed  a  dynamic  equilibrium  of  aggregation  and 
disaggregation,  the  model  in  ^^  requires  that  the  disaggregation  rate  be  modulated  by  the 
amount  of  plaque  via  a  feedback  mechanism  that  minimizes  the  changes  that  occur  in  the 
brain.  While  some  feedback  mechanism  -  for  example,  a  cytotoxic  T-lymphocyte  response 
""^  -  may  exist,  our  model  shows  that  when  A/3  production  is  incorporated,  a  steady-state 
A/3  burden  can  be  achieved  by  a  constant  fragmentation  rate  (i.e.,  in  the  absence  of  any 
feedl)ack). 

Nonetheless,  our  estimate  for  r  is  of  the  order  of  10"'  (rather  than  10"^  or  smaller),  and 
so  it  is  conceivable  that  certain  aggressive  forms  of  the  disease  in  humans  or  mice  (strongly 
amyloidogenic  presenilin  mutations)  may  achieve  a  polymerization  ratio  greater  than  1.  More 
importantly,  while  studies  '^"^--^  assess  Af3  burden  by  counting  plaques,  a  more  recent  study 
measures  AP  burden  biochemically,  and  these  more  refined  data  suggest  that  A/3  burden 
is  correlated  with  clinical  severity  of  AD,  which  is  not  inconsistent  with  linear  growth  of 
A3  burden  (i.e.,  the  supercritical  regime  in  our  model).  Although  further  research  may 
be  required  to  further  elucidate  whether  A/3  burden  is  in  a  steady-state  or  is  experiencing 
continual  growth  in  symptomatic  AD  patients,  our  identification  of  two  regimes  (steady- 
state  and  supercritical)  suggests  that  there  are  three  types  of  possible  treatment  outcomes: 
a  reduction  in  A3  burden  from  a  pretreatment  steady-state  to  a  post-treatment  steady-state 
(Figure  1),  a  change  in  regime  from  pretreatment  growth  to  post-treatment  steady-state 
(Figure  4a),  and  a  reduction  in  growth  rate  from  a  pretreatment  supercritical  regime  to  a 
post-treatment  supercritical  regime  (Figure  4b),  These  results  suggest  that  the  failure  of  a 
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drug  to  reduce  the  A/i  burden  in  a  subset  of  mice  or  humans  may  not  necessarily  be  due  to 
drug  inactivity,  but  rather  could  signal  a  pretreatment  supercritical  regime.  Also,  even  if  the 
pretreatment  polymerization  ratio  is  less  than  1,  a  patient  may  revert  to  his  pretreatment  A/3 
burden  within  several  months  after  treatment  is  discontinued  (see  Figure  5b);  consequently, 
it  is  important  to  attempt  to  measure  the  A/3  burden  throughout  the  course  of  treatment  in 
a  cliniial  trial.  In  contrast,  if  the  pretreatment  polymerization  ratio  is  greater  than  1,  then 
after  discontinuation  of  treatment  a  patient's  A/3  burden  will  always  be  smaller  than  if  he 
had  not  received  treatment. 

Equation  (6)  provides  a  simple  formula  for  the  steady-state  A/3  burden  in  terms  of  the 
polymerization  ratio  r  (provided  r  <  1)  and  the  production  rate-to-loss  rate  ratio,  pi  I.  As 
pictured  in  Figure  1,  the  A/3  burden  is  a  highly  nonlinear  (increasing,  convex)  function  of 
the  polymerization  ratio,  and  approaches  infinity  as  r  approaches  1.  Hence,  the  impact  of 
a  modest  reduction  in  r  appears  to  be  sufficient  to  produce  a  many-fold  reduction  in  the 
A/3  burden;  for  the  parameter  values  in  Table  1,  a  40%  reduction  in  the  production  rate 
via  a  7-secretase  inhibitor  achieves  an  18-fold  reduction  in  A/3  burden.  Such  a  reduction  is 
consistent  with  recently  reported  data  in  mice  ■^'  *  and  appears  to  be  sufficient  to  convert 
the  disease  process  to  a  non-pathological  state  2,3,29^  assuming  neuronal  dysfunction  is 
concommitantly  decreased.  Figure  1  also  shows  that,  while  the  magnitude  of  the  A/3  burden 
reduction  increases  with  the  strength  of  treatment,  there  are  decreasing  returns  to  scale  as 
the  treatment  gets  stronger.  Finally,  our  analysis  predicts  that  the  percentage  reduction 
in  A/3  burden  achieved  by  a  given  treatment  is  a  unimodal  function  (i.e.,  first  increasing, 
then  decreasing)  of  the  polymerization  ratio  r.  For  r  <  1,  the  percentage  reduction  in  A/3 
burden  increases  with  r,  and  hence  one  might  expect  that  a  larger  percentage  reduction 
would  be  achieved  in  forms  of  AD  that  are  particularly  aggressive  (e.g.,  the  Swedish  familial 
AD  mutation  ^',  the  trisomy  21  nnitation  that  occurs  in  Down's  syndrome  ''■^).    However, 
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as  explained  in  the  previous  paragraph,  the  A/3  burden  continues  to  grow  in  the  face  of 
treatment  if  the  post-treatment  polymerization  ratio  is  greater  than  1. 

Our  model  is  flexible  enough  to  incorporate  a  variety  of  treatment  approaches  that 
affect  any  of  the  four  model  parameters.  Equation  (5)  allows  us  to  perform  an  apples-to- 
apples  comparison  of  production  inhibitors,  elongation  inhibitors,  clearance  enhancers,  and 
fragmentation  enhancers.  For  a  given  percentage  change  in  parameters  caused  by  a  treat- 
ment, this  analysis  shows  that  the  production  inhibitors  (clearance  enhancers,  respectively) 
are  somewhat  more  effective  in  reducing  the  A/3  burden  than  elongation  inhibitors  (fragmen- 
tation enhancers,  respectively).  Equation  (7)  and  Figure  3  provide  an  explicit  comparison 
between  production  inhibitors  and  fragmentation  enhancers,  and  computational  results  show 
that  a  fragmentation  enhancer  attains  its  post-treatment  steady  state  faster  than  a  produc- 
tion inhibitor.  These  comparisons,  coupled  with  future  data  on  drug  toxicity,  may  be  useful 
in  determining  the  optimal  mix  and  level  of  treatments  that  minimize  the  A/3  burden  subject 
to  toxicity  constraints. 

Figure  2  shows  that  it  takes  several  months  for  treatment  to  reduce  the  A.(5  burden  to  a 
new  post-treatment  steady-state,  and  the  return  to  a  pretreatment  steady-state  A/3  burden 
following  the  discontinuation  of  treatment  is  even  more  gradual.  However,  Figure  5b  shows 
that  the  rate  of  A/3  turnover  is  highly  sensitive  to  the  pretreatment  polymerization  ratio  r, 
suggesting  that  aggressive  forms  of  AD  (having  high  r)  are  likely  to  take  much  longer  to 
get  under  control.  Hence,  our  lack  of  a  precise  estimate  for  the  pretreatment  polymerization 
ratio  prevents  us  from  accurately  predicting  the  rapidity  of  response  to  treatment,  but  in 
some  cases  this  response  may  be  considerably  quicker  than  suggested  in  Figure  2. 

This  general  state  of  affairs  is  similar  to  HIV  infection,  as  revealed  several  years  ago 
in  two  seminal  papers  ^^■■'•*.  The  stability  of  plasma  HIV  levels  over  many  years  lulled 
the  research  community  into  believing  that  HIV  was  a  relatively  static  disease  process. 
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but  a  perturbation  of  this  steady-state  by  a  powerful  protease  inhibitor  revealed  a  highly 
dynamic  equilibrium,  where  a  high  virus  production  rate  was  offset  by  an  equally  high  virus 
chvirance  rate.  Hence,  antiviral  therapy  led  to  a  precipitous  drop  in  plasma  HW  levels,  and 
discontinuation  of  treatment  allowed  for  a  rapid  return  of  HIV  levels  to  the  pretreatment 
steady-state.  To  borrow  the  analogy  often  used  in  the  HIV  field,  the  water  level  (i.e..  the 
A/9  level)  in  the  bath  may  change  very  slowly,  but  the  tap  (production  of  A/i  monomer  from 
APP)  and  drain  {A(5  loss)  may  be  operating  at  deceptively  high  rates.  It  is  worthwhile  for 
AD  researchers  to  keep  in  mind  the  lessons  learned  by  the  HIV  clinical  research  community: 
(i)  if  toxicities  of  various  types  of  agents  are  non-overlapping,  then  drug  combinations  (e.g., 
a  7-secretase  inhibitor  plus  an  A/i  vaccine)  are  likely  to  outperform  monotherapy;  (ii)  watch 
out  for  drug  resistance,  which  is  the  Achilles  heel  of  HIV  treatment  '*'';  and  (iii)  given  the 
dynamic  equilibrium  oi  A/3  -  which  implies  that  AD  patients  may  require  chronic  treatment  - 
and  the  interaction  between  Ap  production  and  Notch  signalling  '■^■^-•^''^  efforts  should  focus 
on  developing  drugs  that  avoid  long-term  toxicity,  which  is  beginning  to  plague  the  HIV 
community  ^^. 

Several  of  the  model  parameters,  particularly  the  production  rate  and  the  fragmenta- 
tion rate,  are  difficult  to  measure  in  vivo.  One  benefit  of  our  analysis  is  that  equations  (5) 
and  (6)  allow  an  indirect  estimation  of  the  production  rate  p  and  the  fragmentation  rate  / 
given  pul)lished  data  on  the  A/9  burden  ^  and  the  fraction  of  Afl  that  is  monomer  ^-  ^^  (or, 
more  generally,  any  quantity,  such  as  mean  polymer  length,  that  can  be  derived  from  the 
steady-state  distribution  of  i-mers).  Equation  (5)  can  also  be  used  to  measure  the  strength 
(i.e.,  percentage  inhibition  or  enhancement)  of  a  treatment  in  vivo,  given  data  on  the  pre- 
treatment and  post-treatment  steady-state  Al3  burdens.  Moreover,  a  transient  analysis  of 
our  model  (not  provided  here)  would  allow  for  indirect  estimation  of  a  third  model  parameter 
from  serial  measurements  of  post-treatment  A(i  burden. 
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Our  study  has  two  important  limitations.  First,  our  model  does  not  currently  differ- 
entiate between  the  various  forms  of  A/3  {A(3iQ  vs.  A/342,  soluble  vs.  insoluble).  However, 
while  we  have  been  careful  to  restrict  the  model  to  A/i  polymerization  per  sc  and  ignore 
protofibril-to-fibril  conversion  and  plaque  growth,  it  is  clear  that  our  model  can,  in  fact,  be 
viewed  as  a  representation  of  plaque  growth  if  deposition  onto  plaques  and  fragmentation  of 
plaques  are  primarily  due  to  monomer  addition  and  break-off,  respectively.  While  monomer 
deposition  onto  plaques  is  believed  to  be  the  main  form  of  aggregation  in  vitro  '^^'  '*'',  it  is 
quite  likely  that  polymers  coalesce  in  vivo  ^^.  Similarly,  fragmentation  of  oligomers  ''^  and 
plaques  may  be  more  complex  than  the  monomer  release  that  is  assumed  here.  Hence,  a 
more  realistic  model  of  plaque  growth  might  use  a  variation  of  Smoluchowski's  equation  that 
allows  i-mers  and  /-mers  to  coalesce  (e.g.,  at  a  rate  proportional  to  'i''j^Ci{t)cj{t),  where 
A;  =  2/3  or  1,  depending  on  the  porosity  of  plaques)  and  allows  fragmentation  of  monomers 
from  the  plaque  surface  (at  rate  fi^'^Ci{t))  or  the  entire  plaque  (at  rate  ftCi{t)).  While  such 
a  model  would  drastically  change  the  distribution  of  A/3  ?;-mers  relative  to  the  current  model 
(skewing  the  distribution  towards  huge  polymers  that  represent  plaques),  it  is  unlikely  to 
significantly  alter  the  basic  nature  of  our  results  pertaining  to  the  A/3  burden.  However,  to 
the  extent  that  the  model  parameters  including  the  effect  of  treatment  -  vary  for  A/34o  vs. 
A/342  or  soluble  vs.  insoluble  A/3,  then  these  model  generalizations  would  be  worth  pursuing 
in  the  future. 

The  second  main  weakness  of  our  study  is  that  our  parameter  values  ar(>  not  necessarily 
indicative  of  human  AD  brains:  the  loss  rate  /  was  derived  from  the  brains  of  mice  ^,  and 
the  elongation  rate  e  was  measured  from  synthetic  A/3  in  a  low  pH  environment  i'-  '■*.  Also, 
we  crudely  equated  the  fraction  of  A/3  that  is  monomer  with  the  fraction  of  A/3  that  is 
soluble.  Our  sensitivity  analysis  (Figures  5a-5b)  suggests  that,  regardless  of  how  inaccurate 
our  parameter  estimates  are,  a  modest  (several-fold)  change^  in  a  parameter  value  caused  by 
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treatiiieiit  will  still  lead  to  a  clinically  significant  (i.e.,  at  last  three-fold)  reduction  in  the 
steady-state  k(5  burden,  but  that  the  time  to  achieve  this  post-treatment  steady-state  (i.e., 
relaxation  time)  is  quite  sensitive  to  our  parameter  values. 

Despite  the  model's  simplicity,  our  analysis  elucidates  the  basic  nat  ure  of  the  kinetics 
of  the  .\d  burden  in  the  face  of  etiologic  treatments,  and  provides  a  framework  -  and  a  novel 
metric,  the  polymerization  ratio  -  within  which  to  interpret  the  laboratory  and  clinical 
results  that  are  likely  to  be  generated  during  the  next  few  years.  As  more  progress  is  made 
on  the  identity  of  the  toxic  moiety  ^'  ^^ ,  the  elucidation  of  the  inflammatory  and  neurotoxicity 
cascade,  and  the  parameter  values  and  mechanisms  related  to  A/9  polymerization  and  plaque 
growth,  our  hope  is  that  models  such  as  this  one  will  be  refined  and  generalized  to  provide 
some  guidance  about  the  optimal  way  to  employ  the  emerging  anti-A/3  drug  arsenal. 
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Figure  1:  Steady-state  A/3  burden  and  clinical  dementia  rating  (CDR)  score  as  a  function 
of  the  polymerization  ratio  r,  with  p  and  /  fixed  at  the  values  in  Table  1.  A  CDR  score 
of  ().()  (0.5.  1.0,  2.0  and  5.0,  respectively)  corresponds  to  no  (questionable,  mild,  moderate, 
and  .severe,  respectively)  dementia.  The  relationship  between  steady-state  A/9  burden  (left 
ordinate)  and  CDR  score  (right  ordinate)  is  based  on  Table  1  of'-,  as  explained  in  the  text. 
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Figure  2:  A/3  burden  and  CDR  score  versus  time.  A  7-secretase  inhibitor,  which  reduces  the 
production  rate  by  40%,  is  administered  from  day  250  to  day  C20  (i.e.,  for  1  year). 
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Figure  3:  Post-treatment  A/i  burden  and  CDR  score  versus  the  percentage  inhibition  of 
production  rate  (via  a  7-secretase  inhibitor)  and  the  percentage  increase  in  fragmentation 
rate  (e.g.,  via  an  A/3  vaccine). 
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Figure  4:  A/i  burden  and  CDR  score  versus  time.  An  A/i  production  inhilntor  is  begun  on 
day  30  in  (a)  and  on  day  245  in  (b),  and  is  assumed  to  continue  indefinitely.  In  (a),  the 
polymerization  ratio  /■  is  greater  than  1  before  treatment,  and  is  less  than  1  after  treatment. 
In  (b),  the  polymerization  ratio  r  is  greater  than  1,  both  before  and  after  the  start  of 
treatment. 
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Figure  5:  (a)  The  pretreatment  steady-state  A^  burden  divid(>(l  In-  the  post-treatment 
steady-state  A/3  burden,  as  a  function  of  the  pretreatment  polymerization  ratio  r.  (b)  The 
relaxation  time  (see  the  text  for  a  definition)  to  a  post-treatment  steady-state  k3  burden, 
as  a  function  of  the  pretreatment  polymerization  ratio  r .  In  botli  figures,  it  is  assumed  that 
a  7-secretase  inhibitor  reduces  the  production  rate  by  40%.  In  thi'  parameter  estimation 
procedure,  r  is  varied  by  changing  the  fraction  of  A  J  that  is  monomer  (i.e.,  C\jc). 
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